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■ Abstract. We calculate the off-center observational relations in a spherically symmetric dust universe that 
is inhomogeneous at small redshifts. In contrast to the usual model, in which the CMBR dipole is interpreted 
as a Dopplcr effect due to peculiar velocity, our model explores an alternative interpretation, in which the 
CMBR dipole is non-Doppler, and the observer is comoving with the mean matter flow. We do not assume 

^ ' a background frame relative to which peculiar velocities are calculated. Our analysis is fully nonlinear and 

I the density contrast is not assumed to be small. We obtain exact expressions for the Hubble and deceleration 

' parameters, and find that both parameters have quadrupole anisotropics, but no dipoles. A simple numerical 

\ procedure for calculating the CMBR dipole anisotropy in our model is presented, and the observed 0.1% dipole 

■ is shown to be reproducible with a reasonable choices of parameters. 
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1. Introduction 



Cosmological models with inhomogeneous spherically symmetric regions have been extensively studied (e.g. 
O \ Wesson 1979; Raine and Thomas 1981; Lynden-Bell et al. 1988; Bertschinger and Juszkiewicz 1988; Bertschinger 

and Dekel 1989; Panek 1992; Meszaros 1994; Moffat and Tatarski 1995; Tomita 1995). With some exceptions 
(e.g. Tomita and Watanabe 1989; Paczynski and Piran 1990; Arnau et al. 1994), the inhomogeneity in the 
matter is generally treated as a perturbation of a homogeneous background, which restricts the analysis to small 
density contrast and small peculiar gravitational potential in a linearised theory. The CMBR dipole is then 
I interpreted as a Doppler effect arising from our peculiar motion relative to the background. 

I Linearised inhomogeneous models with a position-dependent Hubble parameter Hq can be used to explain 

why measurements of Hq (see Riess et al. (1995) and Van den Bergh (1995) for recent results) tend to give 
relatively high values from methods dependent on local observations, and lower values from methods based on 
more distant observations. (The Hubble parameter relevant to large-scale cosmology should be low to avoid 
contradiction between the universal age calculated from iJp and Op, and the limits due to the age of constituents 
of the universe). 

In this paper the inhomogeneity is treated self-consistently, without linearised perturbations, thus allowing 
for nonlinear peculiar gravitational potential and density contrast (while also incorporating the case of small 
density perturbations). We therefore consider cosmological models which may be homogeneous on large scales, 
but which can incorporate arbitrary isotropic inhomogeneity on smaller scales (say about 100 Mpc). In order 
to exploit the possibilities of such models, while simultaneously exploring alternatives to the standard models, 
we take the dipole anisotropy of the CMBR, and any anisotropics in the Hubble and deceleration parameters, 
to arise not from the observer's peculiar motion (relative to some background), but from the observer's off- 
center location, i.e. from the inhomogeneous gravitational field. Galaxies follow unperturbed motion in this 
field, so from their viewpoint there is no peculiar velocity. For comoving observers away from the center of 
symmetry, the CMBR dipole anisotropy will be interpreted as a cosmological effect, rather than the effect of 
local peculiar motion. A similar model of the CMBR dipole was given by Paczynski and Piran (1990), and 
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non -Dopplcr dipolcs arising from large scale perturbations were considered by Bildhauer and Futamase (1991), 
and Langlois and Piran (1995). The CMBR itself cannot distinguish between these two interpretations. Indeed, 
the observer's peculiar velocity is not directly observed, but postulated as a reasonable source of the CMBR 
dipolc. Our model thus investigates an alternative explanation. 

As pointed out by Paczynski and Piran (1990), observations have not decisively established the reality of 
large-scale departures from a mean Hubble flow, so that alternative models of the CMBR dipole can help to 
develop new ways of testing the standard model. It turns out that both Hq and in our model (i.e. a freely 
falling off-center observer) have no dipole anisotropy, whereas the standard model (i.e. peculiar motion towards 
a Great- Attractor-like structure) does predict dipoles in Ho and qo {cf. Nakao et al. 1995). Thus, in principle, 
measurements of Hq via the distance rcdshift observations could distinguish between these alternatives. 

Exact formulas for the Hubble and deceleration parameters are given in §4. These arise from the local 
solution of the null-geodesic equation near the off-center observer, as discussed in §2. In §5 we show how to 
calculate the CMBR dipole, and in the appendix the observed value is obtained using a reasonable choice of 
cosmological parameters. 



2. Off— Center Geometry 

The Lemaitre-Tolman-Bondi (LTB) metric (Bondi 1947) gives the geometry corresponding to a spherically 
symmetric space-time with a pressure-free dust matter distribution, prior to the formation of shell-crossing 
singularities: 

ds'' = -de + ^^^^dr'' + R{r, t)\de'' + sin^ 9 d^^) (1) 
1 + / (r) 

where x'^ = t is cosmic proper time and = {r, 9, (p} are comoving. The arbitrary function / determines the 
curvature of the spatial hypersurfaces. The units are chosen such that c = G = 1, and the cosmological constant 
A = 0. Einstein's field equations reduce to 

# = ^+/ (2) 

(3) 



where p is the proper density, and M(r) is the gravitational mass inside the sphere of radius r. A prime denotes 

d/dr, and a dot d/dt. The four-velocity of the dust is u'^ = S^. 

Let Pq label the spacetime event of an off-center observer at time t = to, i.e. the event {t = to,r = ro, 9 = 
7r/2,</. = 0}. 

With fc'' the tangent to the geodesies x'^{v), the geodesic equations k^.^k'^ = take the form 
d^t R'R' fdrV .^^o ^ 

k' dt dr fR" f \fdrV ,^ R ^ 

dH Rd9^^^Rdrd9_ ^^^^^^^ /d^y ^ ^ 
dv^ Rdv dv R dv dv \dv J 

d'^cj) ^Rd(f>dt ^R' dr d<j) 2 cot 6» '^'^ - (7) 
dv"^ Rdv dv R dv dv dv dv 

2 



where L? = {d9/dv)'^ + sin^ 9{d(j)/dv)^. The first integral 



dtV fdrV „2r2 ^ 



^dv J 1 + / \dv^ 

specifies that the geodesies are null. If the choice u = at Pq is made, then the past nuU-geodesics through Pq 
are locally expressible as 

a;«(t;) = ^(x«)„t;" +0{v') (9) 

n=0 
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and the functions R and / may similarly be assumed Taylor-expandable. Throughout the paper a subscript '0' 
denotes evaluation at Pq. The geodesic equations (§-0) may be solved locally by equating coefficients of powers 
of V. The solution is of the form {r2, ra, • • • , t2,t3, ■ ■ ■ , 92,0^, ■ ■ ■ , (/'2, 03 • ■ ■} with each coefficient a function of 
the observational angles and ip a.t Pq. 

Since afhne parameters remain affine under linear transformations, v may be fixed by choosing 



h = -I 



(10) 



This choice forces v to increase down geodesies of incoming light, and is equivalent to setting u^k^^ — 1 at Pq 
(Ellis et al. 1985). Evaluation of equation (p|) to lowest order in v results in 



(11) 



The factor e = (or 1) is determined by the motion of rays at Pq toward (or away) from the center of symmetry. 

Equations (^ and ( pi] ) are equivalent, since geodesies that are initially null must remain null. The four 
components of at Pq may be interpreted as follows: ti fixes the gauge of the affine parameter by equation 
(|l0|), ri specifies the geodesies are null by equation (|ll|), and 9i and 0i determine the 2~sphere of observer 
angles {-i^, at Pq. 

Calculation of d and ip in terms of 0i, 0i, and e, is as follows. Let J^, and be unit vectors in the 
surface {t = to], with pointing radially outward from the center of symmetry, K^^ pointing in the direction 
associated with —d/dO, and pointing in the direction of an observed light ray: 



Vi+7 



R' 



5'i 



= 0, 



V(i + /) [i-R^{ei + <j^l)] 

R' 



1 



R 



Si; 



Let 1? and ip be defined (relative to 0, 0, and e) via 

e = = (j}i^ip^O, 
Then (p is given by cosi^s = /^J^ with 6*1 0, which yields 



9i = ^ i9 = 



Also cosi? — J^Kf^i implies 
Equations ^ and (0) invert to 



(p = ei: + (—1)' arcsin {Ripi)^ 
i9 = arccos {—ROi)^ 



sm ip 



R 



cost? 



R 



(12) 
(13) 
(14) 

(15) 

(16) 
(17) 

(18) 



Now, by spherical symmetry, just one angle is required to describe the off-center observations: the angle ip 
between /'^ and the direction of observation at Pq. Observations with common ip must be indistinguishable. 
Attention may therefore be restricted to the plane {9{v) = Tr/2 = d}, for which tp = ip. 



3. Matching Constraints 

We require that the inhomogeneous models become homogeneous on a sufficiently large distance scale. For 
simplicity it is supposed that the transition occurs on a well-defined comoving spherical boundary without any 
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surface density layer (Israel 1966). The coordinates may be extended beyond the boundary (at r = Xi say) 



and into the external FLRW region. For r < x the inhomogeneous LTB solution is 

/ = 0: i?= (9M/2)^ (t-/3)i (19) 

/>0: i? = A/(cosh?7- 1)//, smhr]-r] = fi {t- (3)/M (20) 

0>/>-l: = M (cos 77-1)//, ri-sinr]=\f\i {t~ (3) /M (21) 

in the parabolic, hyperbolic, and elliptic cases respectively. The function (3{r) gives the big-bang hypersurface. 
For r > X the homogeneous LTB (FLRW) solution is 

/ = 0: i?=(9AftV2)^ (22) 

/ = a(2M)^>0: R {cosh -q ~ 1) {M/4a^)^ , sinhr] - i] = 2ah (23) 

> / = -a(2Af)^ > -1 : R = {1 - cos rj) (M/Aa^)^ , 77 - sinry = 2a^t (24) 

where a is a constant, and the large-scale FLRW Hubble and density parameters are given by 

/ = : n = 2/3i, n = 1 (25) 

/>0: H 2ai sinh?7/ (cosh?7 - 1)^ , = 2 (coshr/ - 1) / sinh^ 77 (26) 

0>/>-l: H = 2aism7]/ {l-cosr]f , 17 = 2 (1 - cosry) / sin^ (27) 

The boundary hypersurface has intrinsic metric K^^^, = g^^^, — n^j,n„ and extrinsic curvature n^-i^Kj^Kp, where 

is the unit normal. The Darmois junction conditions (for the case with no surface layer at r — x) &re 

(Israel 1966) that the intrinsic metric and the extrinsic curvature, evaluated on the two sides of r = x, are 
equal. Restricting attention to solutions which have the same characteristic curvatures (parabolic, hyperbolic, 
or elliptic) on both sides of r = %, the conditions are 

aU / : lim [M{x + S) ~ M{x ~ S)] = (28) 

f = : /3(x)=0 (29) 

/>0 : /3(x)-0, /(x)-«[2M(x)]^ (30) 

0>/>-l : /3(x)=0, /(x) = -a[2Af(x)]^ (31) 



Papapetrou (1978) successfully matched an interior collapsing (/ < 0) LTB solution to the Einstein-de Sitter 
model (for which / = 0). It may also be possible to obtain an internal LTB solution consisting of several 
zones with different solution branches (/ < 0, / = 0, / > 0), obeying the Darmois conditions on their concentric 
spherical interfaces. These are however highly non-trivial models, in which the avoidance of shell-crossing 
caustics is a difficult problem. (When dust shells cross the zero-pressure condition of LTB is violated). An 
analysis of some of the issues relating to the matching of LTB solutions will be the subject of subsequent 
research. 



4. Off— Center Observations: Galaxies 

In this paper, idealised cosmological observations are assumed. In particular, it is supposed that area distances 
(or equivalently luminosity distances) and number counts are known as functions of the redshift. Hence any 
statistical, averaging, smoothing, and distance-dependent selection effects are neglected. As is well known 
(Maartens et al. 1996) LTB models are completely specified by complete measurements of the area distance 
and number count as functions of the redshift. In what follows, these two functions are expanded about Pqj up 
to the first two non- vanishing derivatives in the redshift. 



The redshift z as measured at Pq is given by 1 + 2; = u'^fc^ (Ellis et al. 1985), and by equation (10): 

. = -^-1 (32) 
dv 
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The observer area distance D at Pq may be obtained by considering an observational coordinate system (Ellis 
et al. 1985) at Pq. The observational coordinates are = {w, v, i}, cp}, where {w = constant } are the past light 
cones from the observer's world-line at Pq. The metric tensor in observational coordinates is then 

Thus the angular part of the transformed metric g§^,Ch£<py g-ffv on the past-light cone of Pq can be calculated 
locally from the null-geodesic solution and equations (18). Then D follows from the formula (Ellis et al. 1985) 

sin^ '& = gmg^^ - {g^^f' (34) 
Equation (^^ and the local null-geodesic solution lead to 

'dD\ 1 



— (35) 

y (, C1+C2 cos2 Tp ^ ' 



(P-D\ C3 +C4I cosV-l +C5COS2 V +^61 cos-^P 



J 9 / ^ (36) 

rf^Vo (Ci+C2COs2l^)' 

where the constants C/ are given in the appendix. 

The Hubble parameter of a spherically symmetric inhomogeneous universe is often given two definitions 
(e.g. Moffat and Tatarski 1995): 

Rr = R/R, = R'/R' (37) 

which are to be interpreted as radial and azimuthal expansion rates (generally the expansion is anisotropic since 
the shear is non-vanishing) . The two definitions agree at the center of symmetry, as follows from equation (^ 
and the behaviours of /, M, and R near the center (Bondi 1947): 

i?(r,i) =i?'(0.t)r + O(r2), f{r)^]-f"{0y+O{r'^), M{r) ^Im'" {0^ + 0{r^) (38) 
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Using equations (37), we obtain 

(0, t) - i?'(0, t)-^^M'"iO)/6R'iO, t) + i/"(0) = Hr (0, t) (39) 

Moffat and Tatarski (1995) have posed the question of which one is the analogue of the FLRW Hq, when 
off-center. We resolve this issue by the use of a covariant definition: 



where SI is the proper distance from Pq to neighbouring worldlines. This leads to (Ehlers 1993) 

Hn = 



(41) 



where 



e = u^.„ = — + 2- (42) 
R' R ^ ' 



is the rate of volume expansion, and 

C^t^y = - ^6 {gf,u + U^,Uy) + U(^f,U^).xU^ (43) 

is the rate of shearing. For the LTB metric (eq. [|] ) , this gives 
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with all other cr'^ vanishing. Equations ( ^1 - ^ ), (p|), and (|T8|) result in 

Ho = Ci+ C2 cos^ V (45) 

Definition (^) is equivalent to the (covariant) definition of Hq as the initial slope of the redshift-area distance 
curve: 

dz 



as may be seen from equation (|3^). The deceleration parameter may be covariantly defined as (c/. Partovi 
and Mashhoon 1984; Kantowski et al. 1995; Dabrowski 1995) 



and using equation ( |36| ) this results in 



9o = -^o(^) -3 (47) 



C3 + C4\cOS1p\+C5COS^tlj + Ce\cOS1p\^ 

go — 9 o (4ol 

(Ci+C2Cos2V;)' 

The above formulae for Hq and qo are exact covariant generalisations of the FLRW definitions. Note that neither 
of the expressions possess a dipole. In principle, they provide an observational test of these models against the 
peculiar velocity-type models which do predict dipoles. Typical values of the Cj for Great-Attractor data are 
given in the appendix. 

The galactic number count N (per steradian at Pq) is in general dependent on the direction of observation. 
Furthermore, the total integrated number count over all angles and out to some distance, will depend on the 
distance definition used (redshift, area-distance,- • •). One may obtain N (considered as a function of d and (p) 
fron^ equation (19) of Ellis et al. (1985): 

dN ^ snD'^{l + z)dv (49) 

where n is the proper number density of galaxies and s is the selection factor (mean fraction of galaxies that are 
detectable). In general s varies down the nuU-geodesics; here it is assumed constant since such observational 
problems are not the focus of this paper. The function n is determined by n = p/m, where m is the mean 
galactic mass (assumed constant) and p is given by equation (||). Equation (^ and the local null-geodesic 
solution lead to 

ydz^J^ V4^m/ (Ci+Cacos^v)^ 

d'^N\ _( s \ + C9I cosV'l +Ciocos2?/' + C11I cosV'p +Ci2Cos'^ V 
^47rTO/ (Ci + C2 cos2 V)^ 



dz^ 







where the constants C/ are in the appendix. The first two derivatives of N with respect to z vanish, by equation 
(^). Note that, like Hq and goi these N-z relations do not contain any dipoles. 
The luminosity distance -Dl is (Ellis et al. 1985) 

DL=D{l + zf (52) 

All derivatives of D, N , and Dl with respect to z at Pq niay be calculated in principle. However, high-order 
derivatives are of no practical interest as they are well beyond the bounds of current observational sensitivity. 
Given the assumption 2; <C 1 on the boundary {r = x}, the model is well-characterised by ideal observations 
expanded as far as their first two non- vanishing powers in the redshift. 

^Equation (19) of Ellis et al. (1985) has an error: the RHS should be multipicd by dv/dy. 
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5. Off-Center Observations: CMBR 



In this section we present a numerical procedure for calculating the CMBR anisotropy, and in the appendix 
the dipole is evaluated in a special case. 

To obtain the dipole in the anisotropy, it is supposed that there exists a well-defined last scattering surface 
from which the CMBR is emitted isotropically, so that we exclude the effect on the CMBR of processes at 
decoupling, in order to focus on the purely gravitational effect of the nonlinear matter inhomogeneity (c/. Raine 
and Thomas 1981; Tomita and Watanabe 1989; Arnau et al. 1994). We also ignore the anisotropics arising 
from inhomogeneous regions other than our own (in the 'Swiss cheese' model, see Arnau et al. (1994) for 
detailed calculations). The calculated dipole anisotropy is thus due to the off-center location of the observer 
in the inhomogeneous region (figure 1). We do not calculate the quadrupole (or higher multipole moments), 
as the inhomogeneous-matter contribution to the CMBR quadrupole is expected to be swamped by effects at 
decoupling (see Nakao et al. (1995) for estimates in the case of an observer with peculiar motion). 

The dipole anisotropy may be calculated for an off-center observer by noting that the first two terms of the 
multipole temperature expansion are T and AT cos V'- Calculation of the difference in redshift of CMBR rays 
along {ip = 0} and {?A — tt} then suffices to obtain AT/T. These trajectories are the radial nuU-geodesics: 

m . (53) 

dr Vi+TH 

With (^3|), (|3^), and L = 0, equation (jj) expresses the redshift as 



z{r) 



The coupled first order o.d.e.s, (53) and ( p4[), are readily solved numerically, with initial data {2(^0) — 0, i(ro) — to}. 
In the hyperbolic and elliptic cases (eg. ]20|-|2l|) we do not have R[r,t{r)] explicitly. For numerical purposes, 
we write instead R[r,rj{r)] and solve equations (p3|-p3) coupled with a third o.d.e.: 



= 77 [r, ri{r)] + 77' [r, t{r), 'qir)] (55) 

The cosmological scenario under consideration in this paper is illustrated in figure 1. Let the proper times at 
events A, C, and X be denoted t^, tc, and tx. The redshift between Pq and event X (event C) is zx (zc), 
and Zyix is that between A and X. The FLRW redshift z^c between B and C is given by 

f {tc/tB)^ (/ = o) ] (n^n^) 

l + ^BC=< (cosh77c-l)/(cosh77B -1) (/>0) > = ; , (56) 

[ (cos77c~l)/(cos77B-l) (0>/>-l) J 



and by ts — tA- The parameter t] is defined in equations (|22 

Now the redshift in the ip ^ Q direction, minus the redshift in the opposite direction is 

Az = (1 + zc) (1 + ZBc) - (1 + zx) (1 + ZAx) (57) 
and the exact expression for the temperature dipole anisotropy is 

^ = ^ {(1 + zc)-' (1 + ZBc)-' - (1 + zx)-' (1 + ZAx)-'] (58) 

which is approximately — Az/2 for small redshifts. 

For numerical calculations, the gauge freedom in the radial coordinate should be removed. Assuming a finite 
and non-zero proper density, equation (|^) implies M{r) = 0[r^) for small r. Therefore the coordinate freedom 
in r may be removed by choosing a mass-coordinate: 

2M{r) = (59) 

such that the remaining (non-trivial) arbitrariness in the dust solution (given by two arbitrary functions; Bondi 
1947), corresponds to physically distinct models. 

Equations (p3|-^) give a simple procedure for CMBR dipole calculation. A numerical example is given in 
the appendix, reproducing the observed 0.1% dipole for a reasonable choice of parameters. 
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6. Discussion 



In §4 we obtained that the angular dependences of Hq, qo, cfiN/dz^, and d'^N/dz'^ are functions of jcosV'l, 
and hence have no dipole components in their multipole expansions. Since ideahsed observations taken further 
than their second non-vanishing redshift derivatives {d'^ D / dz'^ , etc.) are well beyond the sensitivity of current 
observational data, it can be concluded that all current idealised galactic observations should reveal no dipole 
dependences, if our model is a valid description of the real universe. If dipoles were observed, our model would 
fall away, and the peculiar-velocity models would be strengthened. 

In the appendix, the main results of this paper are illustrated by application to the Great-Attractor data. 
The Hubble parameter result 

Hq w (45 + 12 cos^ V) kms^^Mpc"^ (60) 

shows a quadrupole term with an amplitude of w 12% of the magnitude of the monopole. We emphasise that 
the existence of a quadrupole component in Hq (and in qq) has nothing in common with 'correcting' Hq for 
any 'peculiar velocity' (such as that defined by the CMBR dipole). Furthermore, the quadrupole cannot be 
removed by any such Lorentz transformation at Pq, since the corresponding correction to Hq is dipolar and not 
quadrupolar (Nakao et al. 1995). The GA data, whilst having considerable standard deviation, do represent 
a reasonable local geometry, and suggest that the angular variations of Hq and qq are not negligible in many 
physically acceptable cosmologies. 

The set of physically independent arbitrary parameters at Pq is {R, R, R' , R" , f, f , f"}. The mass out to 
Pq is given by R, but higher time derivatives are not independent, due to the field equation (H). This set of 
seven parameters can be related to a set of seven independent observed parameters by an algebraic inversion. 
For brevity, the result of the calculation is not given in this paper. Note however that one possible choice of 
seven independent observed parameters is {Ci,C2,C3,C4,C7,Cg,Cii}. 

In LTB models, the quadrupole term in Hq vanishes along off-center comoving world-lines when C2 = 0. 
This reduces to the pair of conditions 

(/3')o = 0, [/'-/(4/M)*]^ = (61) 

which imply (in one sense) that the universe is locally FLRW: the first implies that in the comoving frame, 
the big-bang was locally simultaneous at r = tq; the second implies that either / = (corresponding to the 
Einstein-de Sitter model), or that at r = rg the radial gradients of the gravitational mass M and total energy 
W (W = VTT7, Bondi 1947) are related by 

3W dM 

-dW = — (62) 



-I M 

which is satisfied for all r in FLRW models. Hence if the quadrupole vanishes everywhere, the universe must 
be the FLRW model. 

Wesson (1979) considered an exact self-similar dust solution (a special case of LTB), and obtained expres- 
sions for the observables. This paper may be considered a generalisation of Wesson (1979), since for the special 
case he considered, 

_ 2 

3 



R{r, t) 



fir) = (63) 



his results follow from those here. Note however that here the full angular dependences of the observa- 
tions are given (Wesson considered only the azimuthal and radial directions), and that here the luminosity 
distance-redshift relation and number count-redshift relation are considered, as opposed to Wesson's (equiva- 
lent) magnitude-redshift and number density-redshift relations. 

While the off-center position of the observer removes many of the simplifications imposed by isotropy, 
spherical symmetry remains a special case of inhomogeneity. However, assuming that local inhomogeneities 
are spherically symmetric is a first step towards the study of realistic inhomogeneities and their impact on 
observations. 
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Appendix 



a. Specific Models using Great— Attractor Data 

The Great- Attractor (GA) is generally assumed to cause a perturbation of the homogeneous background which 
lies in the linear regime (cf. Panek 1992). Data for the attractor is typically in the form of a peculiar velocity 
field V. In particular, Lynden-Bell et al. (1988) give 

V = VoSo{S-^)t=to (64) 

where S is radial proper distance from the GA, Sq w 87 Mpc, Vq w —570 kms~^, and we take TCq = 
50 kms~^Mpc~^. Quantitative density estimates are far more speculative. 

In the small density contrast and small peculiar gravitational potential limit of our model, we may identify 

V = (r- Hr) (65) 

V / t=to 



and R ^ S. Note that there is no unique generalisation of equation ( |65| ) to the fully nonlinear case. In our 
model, V is not a peculiar velocity, but is used to fix numerical values in the model. 

Hence we require an inhomogeneous LTB model which is specified by equations (|6^) and (^), and an 
assumed density distribution. Note however, that equation (^^ will be simulated only approximately, since the 
matching conditions require V = a,t r = x- Also equation ( p^ ) is singular at the center, whereas we shall 
require V = there, on physical grounds. One possible inhomogeneous model that fits these observational 
constraints is the elliptic solution given by 

0> fir) =-ar2{l-Ai [exp(-A2rVx') -exp(-A2)]} (66) 
/3(r) = A3 [exp (-A4rVx') - exp (-A4)] (67) 



which matches to the closed FLRW solution ( |24|) at r = X: with r the mass-coordinate [eq. (p9[)]. The constants 
A/ are given by 

Al = 80, A2 = 0.43, A3 ~ 1.96 x 10"s, A4 = 0.25 



and we choose a « 1.32 x 10~^'^s~3, so that flo = 1.1. The proper radius out to r = x is chosen to be 
« 96.9 Mpc. The ratio of the inhomogeneous proper density to the homogeneous background density p/p, is 
plotted in figure 2, and figure 3 gives the CMBR dipole. From the plots, the observed 0.1% dipole occurs at a 
radius of k, 80 Mpc, where the proper density is « 4 x 10~'^°gcm~^. 

The (mathematically simpler) parabolic case is inappropriate for description of the GA (in our model), as 
it predicts a cold spot in the direction of the GA. The problem is that the rate of expansion near the GA is 
greater than in the FLRW region. Interpreted via the peculiar velocity function V ^ this means V is positive 
(in the direction away from the GA). We can see this quantitatively as follows. The density ratio is 

- = 7 7 (69) 

-p (t _ /5 - f r/3') (< - /?) 

We must have (3'{r) < for r < x, to avoid singular densities for t — (3 > Q. From (]29|), /3(x) = 0, hence the 
function (3 must be everywhere positive, and the peculiar velocity 

^ = Wor (^^y (io-/3)~^/3 (70) 

must also be positive. 

It is however possible to use the / = assumption to estimate the observational parameters (C/), as they 
are affected only by the local metric, and there is no evidence to rule out / ~ locally to Pq. Hence we take 
S = R + £, where e corrects for the deviation of S from R near the GA, and is approximately constant on the 
time slice {t — to} near Pq- The total mass out to distance from the GA is w 5 x 10'^'' kg (Panek 1992). 
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This fixes the mass-coordinate at Pq, and the observed parameters C/ may be obtained from appendix (b) and 
by equating spatial derivatives of equations (^) and (^5|) at Pq. The results are 



Ci = +45 
C2 = +12 



knis ^Mpc 





C3 = 


-0.02 


1 


Ci = 


-0.4 


1 


C5 = 


-0.01 




C6 = 


+0.6 



kms "'^Mpc 



(71) 



C7 



+4000 kgkm"^Mpc"\ 



C9 = 

Cio 
Cii 
C12 



-200 

-9000 

-200 

+10^ 

+5 



> kgMpc km s 



(72) 



b. CoefRcients 

The coefficients in equations 

Ci 
C2 
C3 
Ci 

C5 
Ce 

C7 

Co 



ij) are (implicit evaluation at Pq) 




i?2 R2 U.R' RR 
97?4 9i?2/ 24i?'i?3 

4/7? 6/i?' 
6/'i?' &R^R' 



12i?2/' 6i?'i?i? 67?2ij' 6i?27?'2 if'R'R 6fR 6f'R 




R'R^ 
47?3 

'W R'^R"^ 



R^R' R^R' R^R''^ 
f'R lORR'^ 2R^R'R" 



R'^R^ i?3 R'R^ 



R^R'^ 



Cin — — 



&f'R 2R^R" 
R^R'^ ' R'^R^ ' R^R' R^ ^ R^R'^ 

6f'R' 15/'i?'2 12RfR' 18/7?'2 



6fR' GR'^R 30RR'^ GRR'R' 6R^f 



R'^R R'^R 



R^R' 



R^R'- 



r:^r' 



m^iV^ ef'R'R 30R'R^ 9R^f' UR'RR 




R^R'^ R'R^ 
Afk AfR 



f'R' 
RR'^ ^ 

_ 3R^ 
~ i?4 



R^R' i?4 
R"f'R 2RR'R" 



R^R' 
2fR'R" 
^ R'^R^ 

AR^R" 



R^R' 
2fR" 
~ R^ 

6f'R 



RR'^ 
12R^R' 
" R^R' " 
2RR'^R" 



R^R'^ 
i?4 



RR'^ 

6RR'^ 
RR'^ ' 



R^R'^ 
3f'R'^ 



R'R^ 
m^R'^ 



4i?3 
~R^' 

3R^f' 



RR'^ 
f'R'R" 



RR'^ RR'^ i?4 
6/i? 6/'i? 6i?4 
~W ^ Wr' ^ 'W 
f'R 6fR' 2R'^ 

" WR^ ^ i?2i?'2 ^ i?i?'3 " 

107?i?/2 AR^R'R" 



+ 



RR"^ 
3fR'^ 



R^R'^ 
6RfR' 



R'^R R'^R"^ R^R' R'^R'^ R^R' 



R'^R^ 

QR'Rf 
'WW 



2fR" 
' RR'^ ' 

8R^R' \ 
R^R' ] 
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Figure Captions 



Figure 1. Schematic for tlie dipole anisotropy. The frequencies of the two CMBR rays are equal at A and B, 
but are then redshifted differently on arrival at Pq. The dotted lines denote constant-i hypersurfaces. 



Figure 2. Plot of the density ratio p/p against radial proper distance S for the elliptic GA model (eq. |66|-|67|) 
at time t — t^. 



Figure 3. Plot of the CMBR temperature dipole AT/T against radial proper distance S for the elliptic GA 
model (eq.M-ra) at time t — to. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/astro-ph/9602033vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/astro-ph/9602033vl 



